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Abstract 

We develop the Keldysh formalism for the polarization dynamics of 
an open spin system. We apply it to the swapping between two qubit 
states in a model describing an NMR cross-polarization experiment. The 
environment is a set of interacting spins. For fast fluctuations in the 
environment, the analytical solution shows effects missed by the secular 
approximation of the Quantum Master Equation for the density matrix: 
a frequency decrease depending on the system-environment escape rate 
and the quantum quadratic short time behavior. Considering full mem- 
ory of the bath correlations yields a progressive change of the swapping 
frequency. 

1 Introduction 

The characterization and control of spin dynamics in open and closed spin sys- 
tems of intermediate size remain a problem of great interest J| ■ Recently, such 
systems have become increasingly important in the emerging field of quantum 
information processing [2] . The quantum interferences of these systems become 
damped by the lack of isolation from the environment and one visualizes this 
phenomenon as dccohcrcncc. Indeed, the inclusion of the degrees of freedom of 
the environment may easily become an unsolvable problem and requires approx- 
imations not fully quantified. This motivates a revival of interest on previous 
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studies in various fields such as Nuclear Magnetic Resonance (3j , quantum trans- 
port 0] and the quantum-classical correspondence problem [SJ E] with a view 
on their application to emergent fields like the quantum computation [7J |SJ [5] 
and molecular electronics ^3 ^2 El E3 ■ 

The most standard framework adopted to describe the system-environment 
interaction is the use of the Quantum Master Equation, derived from the Liouville- 
von Newman equation |14M15| in a fast fluctuation approximation. Interactions 
with the environment occur at a rate given by the Fermi Golden Rule (FGR) pro- 
viding a dissipative mechanism that could induce a non unitary dynamics into 
the system. An overall (conservation) balance condition is obtained by impos- 
ing a convergence into the thermal equilibrium state. While sufficient for most 
traditional applications, this approximation leaves aside important memory ef- 
fects and interferences in the time domain produced by the coherent interaction 
between the system and the bath which are becoming of increasing interest |lfi| . 

The present work focuses on two spin correlation functions in small open 
systems with environmental interactions under conditions where the dynamical 
feedback effects, that go beyond the Fermi Golden Rule, become relevant. For 
this, we will resort to the Keldysh non-equilibrium formalism which leads to 
an integral solution of the Schrodingcr equation. While this novel situation is 
presented here for the first time, the formalism already inspired original ex- 
perimental and theoretical developments in coherent spin dynamics involving 
quantum interferences in the time domain. In particular, it was used to develop 
the notion of polarization waves leading to mesoscopic echoes ^| El > to es- 
tablish the influence of chaos on time reversal (Loschmidt echoes) ^3 120] and 
to establish the possibility of a spin projection chromatography [21]. A rough 
account of many-body decoherence enabled the interpretation of anomalies in 
spin "diffusion" experiments as a manifestation of the quantum Zeno effect \22\ . 
We now make a leap forward in the development of this formalism by showing 
how it deals with open systems. The application to a case with an exact ana- 
lytical solution [23 and where more standard approximations can be obtained 
[3] will show the potential of our proposal. 



2 The Keldysh formalism for open systems 

In this section we make a brief introduction to the Keldysh formalism, summa- 
rizing our results for closed systems of Ref . [21] • Our aim is to extend them to 
open systems. Let us start considering a system with M spins 1/2. The spin 
correlation function, 

m '" u (^1^(0)^(0)1^)' 

gives the amount of the z component of the local polarization at time t on ra- 
th site, provided that the system was, at time t = 0, in its equilibrium state 
with a spin "up" added at the rc-th site. Here, S^(t) = e im/h S^e- im/h is 
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the spin operator in the Hcisenberg representation and \^ eq ) = 12n a N ^eq^ 
is the thermodynamical many-body equilibrium state constructed by adding 
states with different number N of spins up with appropriate statistical weights 
and random phases. The Jordan- Wigner transformation (JWT) S+ = 

exp{i7r X^mJi c^c m }, establishes the relation between spin and fermion oper- 
ators at site n. Symbols and c n , stand for the fermion creation and destruction 
operators, and S„ are the rising and lowering spin operator = S% ± iS%, 
where S 1 ," (u = x, y, z) represents the Cartesian spin operator. The initial 
polarized state is described by the non-equilibrium state l^n.c.) = c+ l^eq.) 
formed by creating an excitation in the n-th site at t — 0. Its further evolution 
is contained in the particle density function j5HJ |23 in the Keldysh formalism 
^mn(*2,ii) = w (*n.o.| c^(^i)c„(t2) l^n.e.) which can be split into contributions 
Gm i : n } (*2, t\) from each subspace with N particles (or equivalently TV spins up). 
Considering that we are in the high temperature regime, i.e., ksT is much 
higher than any energy scale of the system, this enables us to re-write Eq. 

as mi 



P m , n {t) = fG< tm (t,t)-l, with (2) 

N=l 

Notice that the non-equilibrium density G^ m (t,t) depends implicitly on the 
index n that indicates the site of the initial (t = 0) excitation. The expression 
for this initial condition is 

<r (0, 0) = { (§EH,< + M • (4) 

Here the first term is the equilibrium density and it can be seen that is identical 
for all the sites. The second term represents the non-equilibrium contribution 
where only the n-th site is different from zero. In general, this density function 
evolves under the Schrodinger equation expressed in the Danielewicz form |28j . 
which becomes: 



l,k 

+ E r r ^ ) (*2,t*)E fc 5 w (t Jfe ,t J )G^ r) (*«.*i)d**<tt«. (5) 

l,k J ° J ° 

are the exact retarded (t 2 > t\ > 0) and advanced 
(0 < £2 < ti) two particle propagators or Green's functions of the many-body 
system. 
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The first term in the rhs of Eq. JSJ can be seen as a generalization of 
the integral form of the (reduced) density matrix (p(t) — e _lWt / R p(0)e lW */ ? ') 
projected over a basis of single particle excitations in its real space represen- 
tation. This term is all one needs to solve systems such as a finite or infinite 
one dimensional chain with nearest neighbors XY interaction In contrast, 
systems with topological defects [23], long range interaction or Ising terms in 
the spin Hamiltonian present complex many-body effects in the particle descrip- 
tion. These lead to mean-life, ImE R , of the single particle states, producing 
the non-conservation of probability on the retarded and advanced propagators, 
G R and G A . In this case, the second term would collect incoherent reinjections, 
given by £<, that compensate any eventual "leak" from the coherent evolution. 
They also can account for processes not conserving the spin projection. A key 
idea in this paper is that a similar effect of density non-conservation appears 
when one attempts to reduce the whole XY system into a "system" of 2 spins 
and an "environment" with M — 2 spins. Under these conditions the sum in Eq. 
||3J| will run only over the subspaces allowed within the "system", N = 1,2. The 
effects of the "environment" will be included in the form of self-energy terms, 
^k~l(tk,ti) and T, R i(tk,ti) modifying the reduced "system". 

If we replace Eq. (J5J into Eq. © and perform the summation in the ./V 
index only over the "system", the result can be seen as the sum of two con- 
tributions reproducing the structure of Eq. (JSJ. Then, the first term will be 
called the coherent contribution because it is related to the initial condition 
within the "system" . The evolution of this initial density decays with time t as 
a consequence of its escape towards the region called the "environment" . 

The second term will account for the thermodynamical nature of the "envi- 
ronment" when M — > oo. It can be seen as a boundary condition that modifies 
the density of the "system" . If the mean occupation of the "environment" is 
lower than that of the "system" , there will be a flux of probability from the 
"system" to the "environment" included in the formalism through the retarded 
and advanced propagators. On the other hand, if the "environment" mean oc- 
cupation were higher than that of the "system" , it would establish a probability 
flow from the "environment" to the "system" and this could be seen as if the 
"environment" were injecting probability into the "system". The evolution of 
this injected probability is described by the second term in Eq. © which will be 
called the incoherent contribution. Thus, the probability within the "system" 
would be fed by the "environment" . 

In general terms Eq. ((2J and © involve two time functions. In order to 
take profit of the information hidden in the time correlations, it is convenient to 
use the new time-energy variables [i, e] • This is inspired in the Wigner co- 
ordinates that exploit the spatial correlations to define the position-momentum 
variables [a;,^]. Appendix A shows how this procedure is performed. Applying 
this technique to Eq.© we obtain 
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^m,m V ' J — 

/oo poo A A 

/ E GS , (^+¥)</ N) (0-0)G^ , (e-f)expH^)^^ 
-oo J-oo k,l=-l 

ft poc poo _____ 

o J —oo J —oo ^ 

exp{-icj(t - ti)> — — d*;. (6) 

We will apply this formalism to a system of M spins I arranged in a chain. 
Their interaction through an XY coupling enables the swapping between nearest 
neighbor spins. In the presence of a magnetic field, the Hamiltonian is 

M M-l 

H chaln = 2_ Su^n ~ \ +5 __] ^n,n+l (&n ^n+l + C.C.), (7) 

n— 1 n— 1 

which has a Zeeman part, Tl z , proportional to S* with Ml n the Zeeman energy; 
and a swapping (flip-flop) term, H. XY , where J n n+1 is the coupling between sites 
n and n + 1. 

This simplified model can be used as an approximation to real 13 C— 1 H 
systems in an NMR cross-polarization (CP) experiment |3U][3]. We will model 
the 13 C and 1 H nuclei, close to the Hartmann-Hahn condition, as the first two 
sites of the linear chain and the rest of the chain would represent the proton 
spin bath or "environment" . 

Instead of solving a high dimensional spin Hamiltonian Q describing the 
"system" plus the "spin bath" , the JWT provides a map into a fermionic sys- 
tem. For a one dimensional chain or ring with nearest neighbor interactions the 
dimension of the Hilbert space can be reduced from 2 M to M enabling the calcu- 
lation of different aspects of spin dynamics ^Jj^lEI] and quantum coherences 
Since the interaction is restricted to nearest neighbors, the only non-zero 
coupling terms are proportional to c^c n+1 — S„S~ +1 . Each subspace with (^) 

states of spin projection (Y_^ =1 S^J — N — M/2 is now a subspace with N 

non-interacting fermions. The eigenfunctions of these sub-spaces are 

expressed as single Slater determinants built-up with the single particle wave 
functions ip a of energy e a . Under these conditions Gm l ,n } = G^ n for all N. 



3 A two-spin system connected to a spin bath 

We label the "system" sites with the numbers —1 for the 13 C and for the 
■"■H containing the initial excitation. Thus, we want to obtain an analytical 
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expression for the local polarization at each site of the "system" that, according 
to Eq. J2J), is proportional to the particle density Green's function Gf^t) (for 
t = -l,0). 

In these conditions, the effective (reduced) Hamiltonian is: 

where E-x = Ml_ l and Eq = fifl , are site energies and V = 4 represents the 
swapping strength. For simplicity we will take E-\ = Eq. The self-energy 
renormalizes the site energy of the 0-th site due to the action of the rest of the 
chain 32 . This procedure makes possible to trace out all the degrees of freedom 
of the environment without loss of information. It satisfies 

E*( £ ) = (^) 2 Ef (e), 

where Vq is the system-environment coupling through the 0-th site. In a finite 
chain (e) can be calculated using the recurrence relations 



V c 2 

e-E~^ +1 (e)' 



X? (e) = - 5 C ^ R (9) 



stating with Ej^(e) = 0. Here E c is the center of the energy band of the homoge- 
neous linear chain that is acting as the environment, V c is the nearest neighbor 
hopping within the chain. In this case, E^(e) is the ratio between polynomials 
of degree M — 3 and M — 2 on e. The roots of the denominator are the M — 2 
eigenenergies of the environment. This functional dependence accounts exactly 
for the memory effects in the "environment" and describes a variety of inter- 
ference phenomena such as quantum beats and mesoscopic echoes. In order to 
include the effect of irreversible loss of information and simplify the calculations, 
we let the number of spins conforming the linear chain tends to infinity, that 
is, M — > oo. On that situation, Sf (e) = S^_ 1 (e) = E fl and Eq. © becomes 
a Dyson equation 0|. When the energy e lies within the band of propagating 
excitations, |e — E c \ < 2 \V C \ , the solution is 

E fl (e) = A c (e) - iT c (e) (10) 

For the present problem we will work with the assumption that \Eq — E—i \ , V, Vq <C 
V c . This means that the temporal fluctuations of the environment are faster than 
any characteristic time of the adopted model. This approximation allows us to 
consider that ImS^(£) ~ (y-) 2 T c (E ) = T and the level E becomes broaden 
according to the Fermi Golden Rule. Since typically, Eq ~ E c , the correspond- 
ing shift ReS^(e) ~ (^■) 2 A c (£'o) is a small correction that can be neglected. 

In order to obtain the contribution of the coherent term in Eq. ©, we 
need to compute the Fourier transform of the product of two propagators 
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Figure 1: a) Schematic representation of the spin system at time t = 0. The 
shaded region stands for the thermodynamic equilibrium state at high temper- 
ature and establish a background probability level. The black filling represents 
the excess of probability over the equilibrium state which is responsible for the 
observed dynamics. In b) the same system as in a) after the JWT, that is, under 
the particle point of view. Note that in this situation the background contribu- 
tion is removed and the dynamics is described by the excess of probability AP. 
In c) we represent the complementary problem of the case b). Here the black 
filling stands for the hole that represent the excitation. In this representation 
it is easier to calculate the memory effects in the bath (see text). 

The evaluation of the incoherent contribution requires some explanation 
about the model for the (e, fcj) function. Following Eq. (3.15) in Ref. [2H! 

E< (e,U) = i2ro(£)fi(M;)4,o<W, 

with r (e) = T as was previously presented and f„(e,^) = (| + AP) 6{ti) 
stands for the occupation factor of the spin bath. Initially, all the 1 H nuclei 
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are equally polarized and this represents the initial condition at t = for the 
environment. A schematic representation of this situation can be seen in Fig. 

(a). All the dynamics arises from the excess of probability AP at the 1 H 
sites. It is interesting to note, as will be explained later, that the background 
probability (shaded region) does not contribute to the dynamics of the system, 
neither in Po,o(i) nor in P_i i o(t). Having this in mind, the initial condition, 
in the particle language can be expressed with a normalized occupation factor 
i n (e,ti) = 1 x 6(ti) for the 0-th site ( 1 H nucleus) and the spin bath, while the 
— 1-th site is empty of excitation f_i(e,fj) = x 9{ti), as can be seen in Fig. ^ 
(b). At t = we allow the environment to interact with the system and starts 
injecting probability into the system through the 0-th site. Then we have 

Zfi{e,t i ) = i0(t i )2T5k,o5 O j- (H) 

To evaluate Eq. @ we need Gq (t,t) from Eq. © which is determined by Eq. 
1)1 lfl. Taking into account that for a two spin system the sum in cq. has 
only two terms, N = 1 and N — 2, we obtain 

Po,o(t) = 1 - ^ jr^r + ^ irrr cos(crf), 12 

2cos^(0) 2cos^(0) 

where we have defined a = ^ ^/(l — (^/) 2 ) and = arctan | ^ ( ^ ) 2 - 
The same calculations for G'l 1 _ 1 (t,t) leads to 

p_ ll0 (t) = i - ^t|> _ |?iz|} cosM _ 20) . (13) 

2cos^(0) 2cos^(</>) 

For both, Po,o(t) and P-ifl{t), the —1 term in the rhs of Eq. (0) cancels out 
with the sum of the coherent evolution of the first term in Eq. (0J and the term 
corresponding to the injection in the N = 2 sub-space . This justifies Fig. 0(6). 

Note that for T — > the above expressions tend to the dynamics in two 
isolated sites: 

P ,o(*)=Q + ^co S (2vW), (14) 



1 1 



P_ h0 (t) = ( 5 - § cos(2^/fi) 1 = 1 - \cos(Vt/h)\ 2 . (15) 

It is also interesting to note that the characteristic time for the decay of the 
probability, t-i = jl, is exactly twice that of a single site with the same en- 
vironment. The interpretation of this is that due to the symmetry adopted 
(E = E_i) the particle is half of the time on each site being less affected by 
the interaction with the spin bath. 

Figure [3 shows the behavior of the 13 C polarization. It can be seen that it 
reaches the value of 1 periodically, converging to the equilibrium value of 1 at 
the exponential rate 1/t2- The first maximum occurs at a relatively short time 
compared with T2. This feature is used in the spin swap operation by stopping 
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Figure 2: The 13 C polarization (solid line) is composed by the sum of a coher- 
ent part (dash-dotted line) and an incoherent contribution (dashed line). For 
comparison, the evolution obtained with the secular approximation of Ref. |3] 
(dotted line) is shown denoted as MKBE. Parameters: T/V — 0.25 
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rf irradiation (and hence the interaction) at a maximal transfer. The maxima in 
our curves of P_ x (t) are always equal to one because of the symmetry adopted 
(Eq = -E-i). However, only the first maxima of the coherent component decay- 
ing as exp[— t/ra], i-e. about 0.7 for our choice of parameters, would be useful 
in quantum information processing. The incoherent component of the polariza- 
tion, having no definite phase relation with respect to the original state, bears 
no information on the quantum evolution. This can be observed by NMR inter- 
ferometry as done in Refs. |18U19| . In this case the observed polarization at 13 C 
presents high frequency oscillations consequence of the interference between the 
polarization amplitude that survived at the 13 C and the component returning 
after wandering in the H system. This interference would be diminished if, in 
the last CP, one uses the second maximum. 

Another interesting feature of Eqs. I|12l) and (|13|l is that they have zero slope 
for t = as can be seen in Fig EI This expresses that the quantum nature of 
the problem has not disappeared within the present approximation, in contrast 
with the result obtained by using the secular approximation T V, standard 
in NMR calculations [3] . Performing the same approximation as in Ref . [3] , but 
considering XY coupling between the system and the spin reservoir, we obtain 
for the normalized polarization of the 13 C nucleus 

P-ifl® = 1 - CXP{ 2 " fl ~ CXP{ 2 ~ f} «*(2Vt/K). (16) 

Both Eqs. (|13fl and (|16fl are obtained considering the fast fluctuations ap- 
proximation which leads to r c (e) ^constant. However, comparing Eqs. l|13fl 
and i|16|) it can be seen that the main differences between them are the decrease 
of the swapping frequency and the extra phase that result from the Keldysh 
formalism. The frequency decrease is a natural effect of the damping in an 
harmonic oscillator and hence its meaning is clear. The extra phase provides 
the correct quadratic behavior for short times. Both effects would introduce 
corrections up to a 10% if one attempts an estimation of the dipolar frequency 
(here 2V/h) from the first experimental maximum. However, if the frequency 
is evaluated from the FT of the signal it differs from the dipolar one in a factor 
of (1 — (^y) 2 ) 1 / 2 - This can have important consequences when one attempts to 
perform a quantification of the 13 C— 2 H average distances |33j . 

4 Memory effects of the spin bath 

The 13 C polarization, P_i t o(t), m the Keldysh formalism arises from the coher- 
ent evolution of the initial particle density, for which the environment is a "sink" , 
and an incoherent contribution where the bath acts as a particle "source" . This 
can be compared with the complementary framework. Instead of dealing with a 
"particle" problem let us consider it as a "hole" problem (Fig. ^ (b) and (c) re- 
spectively). On these grounds, at t = 0, all the sites are occupied except for the 
"hole" excitation at the —1-th site. See Fig. ^( c ) where the black color stands 
for the hole excitation. At later times this excitation evolves in the system and 
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Figure 3: It is shown the short time regime for the Polarization function for 
the 13 C nucleous obtained within the Keldysh formalism (solid line) in contrast 
with the the evolution obtained with the seccular approximation (MKBE) of 
Ref. [3] (doted line). Parameters: T/V = 0.25 
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also propagates through the reservoir. The "environment" does not have holes 
to inject back into the "system" but those evolved coherently from the initial 
hole (i.e. S< = 0). Here the environment is a perfect "sink". Thus all the 
dynamics would be coherent, in the sense previously explained. If we add the 
result obtained in this case with that of Eq. we obtain a one for all times 
consequence of the particle-hole symmetry. This is a particularly good test of 
the consistency of the formalism because in each result the "environment" is set 
in a different framework. It also shows that the background polarization does 
not contribute to the dynamics. 

This "hole" picture can help us to get a very interesting insight on the 
dynamics in a case where the memory on the environment becomes relevant. 
Consider, for example, the case V = Vq = V c and Eo = E_\. The finite version 
of this effective Hamiltonian applies to the actual experiments reported in Rcf. 
|84| . In this case, the simplifying approximations of the fast fluctuations regime 
are not justified. However, the exact dynamics of the system can be analytically 
obtained if one considers an infinite chain. This enables the use of Eq. i|l(J|) to 
evaluate the propagator in the first term of Eq. JSJl . The integration gives the 
first Bessel function, hence: 

(17) 

A first observation is that the frequency above is roughly increased by a factor 
of two as compared with that in Eq. i|15|) . Since the maxima of P_i,o are 
zeroth of the Bessel function it is clear that the frequency increases slightly 
with time. These are memory effects of the environment that are dependent on 
the interplay between the spectral density of the bath and that of the system. 

We notice that the memory effect can also appear in other condition for 
the bath. For example, if the proton nuclei have random polarizations and the 
density excitation is at site —1, i.e. in Fig. ^ ( a ) f«( e ) = (5) f°r n = 0,1, .. 
representing the 1 H sites filled up to the shaded region; and the 13 C site with an 
occupation i + AP. In this case the excitation propagates over a background 
level (shaded region) that does not contribute to the dynamics. The schematic 
view of this initial condition is equivalent to that of Fig. ^ ( c ) where now the 
black filling represents a particle excitation. The solution of the polarization 
is the first Bessel function, P_i,_i(i) — \-^r Ji(2tV/fr)\ . Apart from the finite 
size mesoscopic effect, this is precisely the situation observed in Ref. al- 
though without enough resolution for a quantitative comparison. The effect of 
a progressive modification of the swapping frequency is often observed in many 
experimental situations such as CP experiments. Depending on the particular 
system, the swapping frequency can accelerate or slow down. Reported exam- 
ples are Fig. 5 on Ref. an( i Fig- 4 on Ref. 35 . This simple example solved 
so far shows that environmental correlations have fundamental importance in 
the dynamics and deserve further attention. 



P-1,0® = 1 
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5 Conclusions 



Summarizing, we have solved the Schrodinger equation within the Keldysh for- 
malism with a source boundary condition which results in an injection of quan- 
tum waves without definite phase relation with the initial state. The model 
proposed allowed us to consider the effect of the environment over the system 
via the decay of the initial state followed with an incoherent injection. We ob- 
tained analytical expressions for the polarization of each of the components of a 
13 C— H system coupled to a spin bath, improving the result obtained through 
the application of the secular approximation [3] in standard density matrix cal- 
culation. 

Of particular interest is the inclusion of temporal correlations within the spin 
bath in a model which has exact solution. On one side, it enabled us to show a 
novel result: memory effects can produce a progressive change of the swapping 
frequency. On the other side, this results will serve to test approximate methods 
developed to deal with complex correlations. 

In general, our analytical results based in the spin-particle mapping, allow 
a deeper understanding of the polarization dynamics. They may constitute a 
starting point for the study of other problems, such as different topologies [23"| 
with XY interaction and the extension to dipolar and isotropic couplings. 
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Appendix A 

Let us define the function 

G q (w> (^2; £i) — / / Gq ^ {t%, f 1Tl )S m ,^' > (t m ,t n )G n ''Q ) (t n , ti)dt m dt n , 

n,m J ° J 

(18) 

which is the second term in Eq. JSJ • A similar expression holds for the coherent 
part. The manipulation that follows is independent on the subspace index (iV), 
and we will keep it implicit. Rewriting the integrand in Ea. (|18|) as 

^0,m(^2, tm^m.nitm, tn)G n Q (t n , t\) 

= ^O.mfe — trn, 2 \ m )^m,n(^™ ~ 2" ™ )^n,o(^« — ™ 2" 1 ) 

Go. m ( £ Ri t2+ 2 tm )S^ n (e , tm +*™ )G^ (ea, — J* 1 ) 



exp [-ie B (t a - t m )/%] exp [-ie'(t m - t n )/h] exp [-ie A {t n - h)/H] ^fl^fff- 

(19) 
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Let's define the macroscopic time as t = hfa + ti) and the quantum cor- 
relation time St = t% — t\ which have related time scales of the injection pro- 
cesses as ti — i (t m + t n ) and St{ — t m — t n . These time scales are associated 
with s = ^(sr + Ea), the energies characterizing the quantum correlation, 
and lj = i (er — sa) the frequencies in the observables. The argument in the 
exponential function becomes 



£Rt2 - £RXm + SAtn ~ £ Ah 
— ErI + Er^j — £Rti — £r 

= hx>{t — ti) + eSt — eSti, 



Rt + ErJj ~ £ R~ti ~ £ R%~ + £ Ati — £ A^J- ~ £ At + SAY 



and also 

s'(t m —t n ) = e'Sti. 
The Green's functions take the form 

(~<R ( c t 2 +t m \ _ /~tR. I , huj t+tj , st+su \ 
'~'0,m\ t R> 2 ' ^0,171^ t ~2~' 2 ' 4 > 

(~iA ( t n +ti \ _ /~iA (- hu) t+tj _ St+Stj \ 
u n,0\ £ A, 2 ) — u n,0\ £ ~ ' ~^T~ 4 > 

y < (J t m +t„ \ _ y < ( I j.\ 
^m^X 1 2 I ~ m,n\ c > "%)• 

Finally due to the fact that the transformation have the property that its Jaco- 
bian is equal to one, we have dt m dt n = dtidSti and de^de^ = hdeduj. Replacing 
all these expressions in the integral of Eq. i|18|) we have 



rt2 pti 

y ' / / G$ im (t2,tm)^m,n{tni}tn)Gn : o(tn:ti)dt m dt rl 



X 



C R (cJ-Hi£ t+li. j_ St+SU_ "\V < (sr 1 + \C A — t+tj St+Stj 
'-'O.m^T 2 ' 2 ' 4 ) l -'m,n\ £ J l i) Kj nfi\ £ 2 ' 2 4 

cxp{-i[fiu(i - ti) + eSt - eSU + e'5U]/h}dUdSU-^-^-^-. 

Ann Z7m Z7r 

Then, we can express Eq. (|18|) as 

G < inc ' (t, St) = J G<™- (e,t) exp [-ieSt/h] ^ , 
and using the last two expressions we can identify 
^ofo mC ' ( £ > *) = 

/ / / / n R (<? _i_ &w t+tj _i_ st+su w < / / i \f-iA /_ hu t+tj st+stj \ 

I I I I °0,m l £ + ~2~ ' ~2 1 4 >^m,n\ £ J l i) U n,0\ £ 2~' ~ 2 I / 

x exp [— i[(hiu(t — U) — eSti + e'StA /h] dtidSti————. 
A similar expression holds for the coherent term. 
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Integrating in energy we obtain G^ (t,i). If we consider that the system 
Hamiltonian is time independent the last complex expression simplifies to 

/t /"OO POO 
/ / ^ ^Q.m( £ + %0 S m,ri( £ '^)G'K.o( £ ~ %0 
.0 J-ooJ-oo n m 



... do; de 

exp [-jw[t - U)\ -z--^-rdti 
Z7T Ann 



which is similar to the second term in Eq. 
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